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We study the thermodyanamics of the two flavor massless Schwinger model on a torus at a finite 
chemical potential. We show that the physics only depends on the iso-spin chemical potential and 
there are marked deviations from a free fermion theory. We argue that spatial inhomogeneties can 
develop in the system at very low temperatures. 
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I. INTRODUCTION 



A study of QCD in the presence of a finite chemical potential is important for our understanding of quark matter 
at finite density [T] . Lattice QCD provides a non-perturbative approach to this problem and there has been extensive 
work performed on this topic [21 [3] . The fermion determinant in a fixed gauge field background is complex (it can be 
made real by summing over a pair of gauge field and its complex conjugate but the result is not necessarily positive) 
and therefore suffers from the so-called "sign problem" . QCD wth a finite isospin chemical potential does not suffer 
from the "sign problem" and the physics of this model has been explored [H E] . 

The Schwinger model (QED in two dimensions) has played the role of a very useful toy model for QCD in four 
dimensions. Generalized Thirring models have been studied in detail at finite temperature and finite chemical po- 
tential [HI [7] ■ One main result at is applies to the Schwinger model is the independence on the chemical potential. 
This can be seen as a consequence of the integral over toron fields [H |9] in the path integral formalism. The issue at 
hand is imposing Gauss's law in the path integral formalism |13j . Imposing Gauss's law in the Hamiltonian formalism 
results in the condition that the time- like component of the electromagnetic potential vanish at spatial infinity [lOj . 
This amounts to setting the toron field in one direction on the torus to zero in the path integral [TT]. This would 
allow for states with net total charge to be present but would break the U{1) global symmetry associated with the 
Polyakov loop in the time-like direction placing the theory in a deconfined phase. We will study the two flavor massless 
Schwinger model on a finite torus in the presence of a chemical potential. We will integrate over the toron fields. As 
expected, the theory will be independent of the chemical potential that couples to the total charge but will depend 
on the isospin chemical potential. 

We start with a definition of the model on a finite torus and state the result for the fermion determinant in a fixed 
gauge field background using the zeta-function regularization 12 . We will address the integration of the fermion 
determianant over the toron fields and show that the result is independent of the chemical potential that couples to 
the total charge. We will proceed to address the physics of the isospin chemical potential. 



II. THE GRAND CANONICAL PARTITION FUNCTION 



A. Model basics 



Let I be the circumference of the spatial circle and let /3 be the inverse temperature. We will use I to set all scales 
in the theory and define r = ^ as the dimensionless temperature. The physical gauge coupling is set to j where e is 



dimensionless. 

The Hodge decomposition of the U(l) gauge field on a Z x /? torus is 



Ai{xi,X2) = '^^^ + dixixi,X2) - d2(j)ixi,X2) -^^X2 

^2(2:1, 2:2) = '^^^ + d2Xixi,X2) + di<j){xi,X2), (1) 
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where — 5 < /i^^ < 5 are the toron fields and x(xi,a;2) generates gauge transformations. The electric field density is 

lirk 

E{xi,X2)^-jj+d^c^{xi,X2), (2) 

where <^(a;i, x^) is a periodic function on the torus with no zero momentum mode and k is the integer valued topological 
charge. The gauge action is 

o 2 jL2 72 I- 

s, = ^ + 1, y d^x{d^4>f. (3) 

The determinant of a massless Dirac fermion is zero unless fc = and the determinant for fc = using zeta function 
regularization is [71 IT^ 

2 00 

^ ^ ' ni,n2 — — oo 

where qi is the integer valued charge of the fermion and ^^y^ is the chemical potential. The Dcdekind eta function, 
?7(r) is given by 

oo 
n=l 

There is an infinite normalization factor that has been removed by zeta function regularization from the above formula. 
This factor does depend on t but only in a trivial manner as to shift the zero point energy. 
We define the Fourier components of (j){xi, X2) according to 

f 00 

'^(2:1,^2) - . e^(^-^+'=^-^)0(fci,fc2), (6) 

47r T3 ^-^ki,k2= — oc 

with 0(— fci, — fc2) = 0*(fci, fe) and the prime over sum implies that fci = fc2 = is excluded. Then 



2 / o-o 2 

' / dMdW = It ^ \4>{h,k2f (kl + ^kl" 



2e- 
and 



(7) 



27r / 87r3r2^fci,fe2=- 



B. Bosonic and toronic partition functions 

We will consider the two flavor Schwinger model with f/i = 92 = 1- We write the partition function as 

Z(^i,^2,T, e) = Zb[T,e)Zt{ni,mu2,T) (9) 

where the first factor is the bosonic (0) partion function and the second factor is the toronic {hfj) partition function. 
Since the total action (gauge and fermionic contribution) is a quadratic function in (p, the bosonic partition function 

is 



^^(-'^) - n„..._ J ,,2 + j^,2) (,2'+ ^ [fc2 ^ (10) 



Starting from Q and after a little bit of algebra, the toronic partition function can be reduced to 



00 



ry , s _4/ N -TTT (m2-mi+2m3+2/i2-i^^j4^) +(mi-iiiV^) +™2 /i i ^ 

^t(Aii,M2,T) = 77 (t) 2^ / ^ d/i2e - / v . ; 2j ^ (^^^■j 

mi ,m2 ,m3 — — 00 ~ 2 
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Consider the integral, 



Z. k 



Ml + 



E 



26 



-•irT(fe+2m3+2/i2- 



.■ fl+f2 \^ 



(12) 



Viewing, z = h2 — i ^^^^^ , as a complex variable, we see that the integrand is analytic in z and periodic under 

z z + 1. Therefore, the integral is independent of iii^^ii^. We explicitly sec that the partition function does not 
depend on the chemical potential that couples to the total charge. Note that the integrand is positive definite, if we 
set (/Xi + /Z2) = but is not the case for a general (/^i + fi2)- One will encounter a sign problem if one tries to compute 
the integral numerically with (yUi + 112) not equal to zero. The integral is the same for all even k and the same for all 
odd k. We can write the reduced integral as 



Z|(t)= ^ y d/i2e-^^('=+2™^+2''=)'; A; = 0,1, 



Setting the dimensionless isospin chemical potential equal to 



jjLi = 27r(/zi - IJL2), 



we have 



Let 



2 '-^ 



mi,m2=—oo 



00 

m2 — — 00 

00 



m2 — — 00 

The final expression for the toronic partition function is 

2 1 CSO 

Ztifij, r) = v-\r)e^ ^ ^ cos {{2m, + k) m/) 6-^-^^^-+^^' Z^{t). 

fc=0 m\=—oo 



(13) 



(14) 



(15) 



(16) 



(17) 



C. Thermodynamic observables 

The only contribution to the isospin number comes from the toronic partition function and is 

ainZ(/x/,T,e) III 
Nr = r — = --/(/.., r), (18) 



where 



. ^ ^ ELo E^,^-oo (2mi + k) sin ((2mi + fc) ^,) e— (^'"^+'=)^Z|(r) 
^ EUE™,=-ooCOs((2mi + fc)M/)e— (2'"i+'=)=Z|(r) 



The dimensionless energy is 

= Eb{T;e)+Et{T;Ni) (20) 



^(r;iV.,e) = - ^^-^(^/'-'^) 



and the contributions from the bosonic partition function and the toronic partition function are written separately. 
The result from the toronic partition function is 

Et{T;Ni)-Et{0;0) 1.2^.^ ^ ^nr^ , ^,2 .2. . 



n=l 
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where 



Xl=o EZ=-oo COS {{2m, + k) ^j) ^ [e— (2™^+'=)^Z2^(r) 



ELo E™,=-oo cos ((2ml + k) (2"u+fe)^z|(r) 



The result from the bosonic partition function is 



(22) 



/2tt 



tanh 



2t:t 



- 1 



- 47r fci 



fci=i 



Trfci 
tanh 1 



(23) 



D. Free fermions 



In order to understand the results for the two flavor massless Schwinger model, it is useful to recall that the partition 
function for free fermions in one dimension is given by 



21 [ I 

\nZf ^ — dp l3p + \n[l 

7^ Jo L ^ 



(24) 



where Hf is the chemical potential for free fermions which we set to ^ in order to be consistent with the two flavor 
notation in ( 14 ) . The free fermion isospin number is given by 



Nf 



dlnZ 



27r' 



(25) 



After subtracting the zero point energy, the dimensionless energy of free fermions at low temperatures is given by 



Ef{T-Nf)-Ef{Q,0) 



- ^1 + 



(26) 



The first term is the fermi energy that grows quadratically with the isospin number and second term is the leading 
order low temperature correction that is positive and quadratic in the temperature. 

III. RESULTS AND DISCUSSION 

We proceed to compare the results for the two flavor Schwinger model with that for free fermions. The result for 
the isospin number in ( 18 ) is plotted in Fig. [l] for several different values of temperature. The linear behavior in 
(25) expected of two flavors of free fermions in (p5|) is the first term in (18) and this is achieved only in the high 



temperature limit where the contribution from the second term goes to zero. The first term is the nai've contribution 
from two flavors of free fermions. The second term is the result of integrating the effect of boundary conditions over 
all possible choices. 

We can use Fig. [ijto see how the isospin chemical potential depends on the temperature at a fixed isospin number. 
The quasi-periodicity seen in the figure is a consequence of /(/i/+27r, t) = /(/x/, r) in ( 19 ). Furthermore, /(/i/, oo) — 0, 
and 



fJ-ijoo) 
2tt 



(27) 



like for free fermions. On the other hand, 



lim /(pt7,r) 



Hi 

2tt 



if < ^7 < TT 

if TT < Pi < 2TT 



(28) 



Therefore, 



\Ni], 



(29) 
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FIG. 1: Plot oi Ni{fj,i,r) versus fii (c.f.(18l) for several different values of r 



for all non-integer values of Nj and \Ni~\ is the ceiling function. The behavior for non-zero and finite temperatures 
is to interpolate between (27) and (29) as shown in Fig. [2j The behavior is shown for values of Nj in the range 
< Nj < 3 in. steps of 0.1. Plots are color coded to show periodicity of f{ni, t). Since f{2mr, r) = for all values of 
T > and any integer n, we see that integer values of Nj are special and behave like free fermions for all temperatures. 

Since the partition function is independent of (/ii + ^J-2), the net charge is zero. But, we can maintain the system at 
a non-zero isospin number, Nj. Since the system can exchange particles with the reservoir, the expectation value of 
the isospin number need not be an integer. Let us assume we start at high temperature with a fixed Nj assumed to 
take an integer value. Since the chemical potential does not change with temperature for this case and remains the 
unique value for this particular value of Nj, the system will remain homogeneous at all temperatures. Now consider 
non- integer values of Nj. Th system will be homogeneous at high temperatures since the chemical potential is different 
for different values of Nj. As the system is cooled and brought down to zero temperature, different values of Nj can 
coexist as long as the different values all have the same ceiling value, \Nj~\ since they all have the same chemical 
potential at zero temperature (c.f.(29)). The system is bound to form inhomogeneities at zero temperature. 

We now proceed to use (jlsj) and (21) to compute the toronic contribution to the energy as a function of the 



temperature at fixed isospin number. Consider the zero temperature limit in order to extract the fermi energy as a 
function of the iso-spin number. Since f{2m:, 0) and g{2mT, 0) are zero it follows that the fermi energy for integer 
values of isospin are given by the free fermion value. As r goes to zero, g{fj,, r) approaches a non-zero limit as long as 
Nj is not an integer. As a consequence, the fermi energy is given by 



Ef{Ni) = [Ni\^ + {2lNi\ + 1) {Ni - [Nil). 



(30) 



and it linearly interpolates between the free fermion values at integer values of Nj. We were unable to analytically 
obtain an explicit expression for the linear coefficient in (| 21|). We numerically evaluated it and found that the second 
term in (21) contributes — , 



and the last two terms in (21) contribute 



3r 

2tt ■ 



The leading behavior of the toronic 



6 




0.5 



1 

T 



1.5 



FIG. 2: Plot of I-ii{t) versus r (c.f.(18l) for several different values of Ni in steps of 0.1 starting from and ending in 3. 



contribution to the energy at low temperature is 

Et{T;Ni)-Et{0;0) 



Ef{Ni 



271 



(31) 



This is qualitativel y di fferent from the free fermion result where the leading term is quadratic in r. The linear 
coefficient of ^ in (31) gets modified to ^ for the total energy when we include the leading contribution from the 



bosonic partition function in (231 



The higher order corrections in t to the energy from the toronic partition function, 



(32) 



is plotted as a function of r for various values of Ni in Fig. [3] for Ni — 0.1, 0.2, 0.3, 0.4, 0.5. Due to quasi-periodicity, 
Er{T, l-Ni)^ EriT,Ni) for (0 < iV/ < 0.5). In addition, Er{T,Ni + 1) = Er{t,Ni). 
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